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Abstract 

We analyze numerically the critical properties of a two-dimensional discretized ran- 
dom surface with extrinsic curvature embedded in a three-dimensional space. The 
use of the toroidal topology enables us to enforce the non-zero external extension 
without the necessity of defining a boundary and allows us to measure directly the 
string tension. We show that a phase transition from the crumpled phase to the 
smooth phase observed earlier for a spherical topology appears also for a toroidal 
surface for the same finite value of the coupling constant of the extrinsic curvature 
term. The phase transition is characterized by the vanishing of the string tension. 
We discuss the possible non-trivial continuum limit of the theory, when approaching 
the critical point. Numerically we find a value of the critical exponent v to be be- 
tween .38 and .42. The specific heat, related to the extrinsic curvature term seems 
not to diverge (or diverge slower than logarithmically) at the critical point. 
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1 Introduction 



The theory of random walks has provided us with a powerful link between statistical 
mechanics and euclidean field theory. An euclidian field theory can be expanded in 
a series of intersecting random walks and a number of rigorous inequalities can be 
proven using the random walk representation. In addition various aspects of univer- 
sality in regularized field theory can easily be understood from the corresponding 
universality of the random walks. It has further been possible to develop a theory 
of random walks, which allows the description of fermions in a geometrical way. 

The theory of random surfaces ought to provide us with even stronger tools in the 
study of string theories. String theory in its first quantized version is nothing but the 
theory of a specific kind of free random surfaces. This has been substantiated during 
the last couple of years, but many questions remain unanswered. One of the greatest 
puzzles is that the formulas derived by Knizhnik, Polyakov and Zamolodchikov 

using conformal field theory seem to make no sense for strings embedded in 
"physical" dimensions (d > 1). In more general terms we do not know how to 
couple two-dimensional gravity to conformal field theories with central charge c > 1. 
On the other hand it seems to be no problem to formulate statistical theories of 
random surfaces in physical dimensions. The use of these random surface theories 
in the context of strings is hampered by our lack of understanding of scaling and 
universality properties of the random surface theories themselves. The situation 
here is quite different from the random walk case. The present work is an attempt 
to clarify the situation for a particular class of random surface models. Let us briefly 
motivate why this class is interesting. 

The class of random surface theories we have in mind is one of freely intersecting 
surfaces with an action having an area term and an extrinsic curvature term. Such 
a model can be viewed as a simplified model of physical membranes. The statistical 
aspect enters due to thermal fluctuations and the statistical fluctuations may drive 
the system to a point where the effective surface tension vanishes and the membranes 
are dominated by their curvature energy rather than their surface tension. This is 
the case for fluid/fluid interfaces and also for the so-called amphiphilic membranes, 
which are membranes formed when amphiphilic molecules are brought into contact 
with water and form bilayers by orienting their polar parts towards water and their 
oily hydrophobic tails away from the water. Of course such surfaces are not allowed 
to self-intersect and from this point of view our models can only be viewed as toy 
models for real membranes. However, since the seminal work of Helfrich || such toy 
models have received a lot of attention 

The random surface models with extrinsic curvature are also interesting as ef- 
fective theories of strings. In this context they were first suggested by Polyakov Q . 



2 



He had in mind an effective string theory which could be equivalent (at least at 
long distances) to QCD with heavy quarks. However, it is now clear that fermionic 
string theories can give rise to effective bosonic string theories which have extrinsic 
curvature terms, and in this context the surfaces should be allowed to self-intersect. 
In the case of the superstring it is possible to integrate out the world- sheet fermions 
IJ. After this integration two types of terms are produced, which both depend on 
the extrinsic geometry of the world sheet : 

Seff = S bosomc + W k (A^) + T -f d 2 i^~ d {(e^) 2 + {D a ntf} (1.1) 

Here g ind refers to the metric induced by target space, nf , i = I, ... ,d — 2 are 
normals to the surface, a = 1, 2 are the tangents and D a denotes the covariant 
derivative with respect to the connection in the normal bundle : 

4 n)Q = r#d a n$ - np a nf. (1.2) 

t is a Dynkin factor coming from the fermionic representation and finally Wk{A^ n ') 
denotes the Wess-Zumino action : 

W k (A {n) ) = — Tr ( - [ d 2 i A A A+ — [ d 3 x A A A A a) (1.3) 
8 V 2 J 3n Jd J 

where we have used the notation A a = A^ a Mij, M y - being the generators of SO(d — 
2) and D a three-dimensional disc bounded by the world sheet. 

One class of terms depends very explicitly on the extrinsic geometry of the world 
sheet and is minimized by smooth, flat surfaces. If K(£) denotes the extrinsic 
curvature of the world sheet we have 

K(0 = + (1-4) 

where ri(£) and ra(£) are the principal curvatures of the surface. Since 

/ d 2 i^~ d (D a nA 2 = J d 2 £^g— d K 2 (0 (1.5) 

this term clearly favours smooth surfaces. The same is the case for the term involving 
the tangents. 

As for the Wess-Zumino-like term W(A), the effect is not so clear. In euclidean 
space-time the term is purely imaginary. It is interesting to compare with the 
situation in the case of random walks. What we have in mind is a "supersymmetric" 
random walk, where we have introduced world-line fermions. In this case it is again 
possible to integrate out the fermions. The result is a Wess-Zumino-like term on 
the world- line and it is the analogue of the term Wk{A) defined in eq. ( |1.3|) . In the 
random walk case the effect of the Wess-Zumino term can be analyzed in detail (see 
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for a continuum treatment and for a random walk approach) . Heuristically the 
result can be described as follows: The amplitudes of back-tracking random walks 
tend to cancel due to the phase factor coming from the Wess-Zumino term and we 
are effectively left with a class of random walks which are much smoother than the 
generic random walks of a scalar particle. In effect the short distance Hausdorff 
dimension of the random walk is reduced from two to one, and the corresponding 
short distance behaviour of the corresponding propagator changes from being like 
1/k 2 to 1/k. As first noticed by Polyakov and proven in the context of random walks 
in [0, the scaling limit of this fermionic random walk leads to a representation of 
the massive Dirac propagator. In this case the effect of world-line supersymmetry 
and the corresponding Wess-Zumino term is clearly to favour "smooth" curves. It 
is natural to conjecture that the two-dimensional Wess-Zumino term ( p..3|) has a 
similar effect. 

As we see all additional terms in ( |1 . 1| ) act in favour of smoother surfaces as 
compared to the class of surfaces singled out by the standard bosonic term Sb OS onic- 
This opens up the possibility of an intuitive and simple understanding of the tachyon 
problem of the bosonic string and its cure in the case of the superstring. The 
appearance of tachyons in a bosonic random surface theory is somewhat different 
from that based on formal continuum manipulations. By definition there can be no 
tachyons in a theory of random surfaces where each surface is given a positive weight 
and which satisfy the principle of reflection positivity. However, the scaling limit of 
such discretized, regularized theories might be pathological from the point of view of 
string theory. This is precisely what happens for the ordinary bosonic string theory 
where the dimension d of target space is larger than one. One can prove that the 
string tension does not scale |12 for such theories and as a consequence the scaling 
limit is not really that of a "surface" theory, but rather that of a theory of so-called 
branched polymers consisting of a minimal surface (depending on the boundary 
conditions) from which the only allowed fluctuations are thin "branches" which 
carry no area. Maybe somewhat contrary to intuition the entropy of such surfaces is 
large compared to that of "smooth" surfaces and completely determines the scaling 
properties of the "surface" theory. An obvious cure is to put in by hand additional 
terms in the action which suppress the "spikes" and favour a smoother class of 
surfaces and this is precisely what is done by imposing world-sheet supersymmetry 
and integrating out the world-sheet fermions as indicated in 

We know there are no tachyons in the superstring theories. However, we also 
know that in the usual continuum description the absence is due to a delicate can- 
cellation between bosonic and fermionic degrees of freedom. From this point of view 
it is intriguing how a discretization of the action ( |1.1| ) and a subsequent scaling 
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limit is able to capture this. The obvious answer would be: by universality. But 
although it makes no sense to derive an effective action like ( |1 . 1|) , except in the 
special dimensions d = 3,4,6 and 10, where a classical superstring theory can be 
defined, the effective action itself makes sense in any dimension and one can con- 
sider the corresponding discretized theory and possible scaling limits. Appealing to 
analytic results there should be a difference between the theories, depending on the 
dimension of target space. At the moment we have no understanding of the physics 
leading to such a difference. The concept of universality triggers a related question: 
How far are we able to modify a discretized version of the effective action (|1 . 1| ) and 
still stay in the same universality class. Due to the imaginary Wess-Zumino term we 
are clearly discussing a rather unconventional class of theories from the point of view 
of statistical mechanics. Nevertheless it does not necessarily mean that one cannot 
apply the conventional machinery of statistical mechanics. The above mentioned 
example of the fermionic random walk provides documentation for this and shows 
that one can indeed talk about universality classes [|7||. It would be very convenient 
if we were able to drop completely the imaginary Wess-Zumino term and still stay 
in the same universality class. At a first glance it looks unlikely, again appealing 
to the fermionic random walk example. One can analyse the random walk where 
one has both an extrinsic curvature term like ( |1.5| ) and the Wess-Zumino term. The 
result is as follows: if we choose the coefficient of the Wess-Zumino term to be zero 
we will in general be in the universality class of the ordinary random walk (and the 
corresponding scaling limit will be that of a free scalar particle)^. As soon as we 
take the coefficient of the Wess-Zumino term to be non-zero we are driven to the 
fermionic random walk, where the scaling limit leads to a massive Dirac propagator. 
The same could be true for the effective action ( |1.1| ). However, since the world-sheet 
is two-dimensional while the world-line is one-dimensional we are discussing vastly 
different theories from the point of view of statistical mechanics. The structure of 
phase transitions is much richer in two dimensions, and although a one loop calcu- 
lation seems to support the conclusion reached in the random walk case jnj, this is 
a purely perturbative argument, based on assumptions which will not be satisfied 
in case we have a non-perturbative phase transition. A priori it is possible that one 
could drop the Wess-Zumino term and still stay in the same universality class. 

The above discussion is closely linked to another discussion of possible critical 
behaviour of membranes (or random surfaces) with extrinsic curvature terms. It 
is believed that one has a transition from so-called "crumpled" surfaces to smooth 
surfaces in the case of crystalline surfaces. Crystalline surfaces can be viewed as 

2 The only exception is if we take a scaling limit where we scale the bare coupling constant of 
the extrinsic curvature term to infinity. It is possible to do this in such a way that we get what 
can be called a rigid random walk p. 
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membranes where the individual molecules have a fixed connectivity, i.e. their 
neighbours are fixed, contrary to the situation for fluid membranes. If the coef- 
ficient of the extrinsic curvature term is zero such surfaces will be "crumpled" (the 
Hausdorff dimension for the statistical ensemble of surfaces is infinite), at least in 
the idealized situation where they are allowed to intersect freely. However, for a 
finite value of the extrinsic curvature coupling there seems to be a phase transition, 
and for couplings above this critical value the surfaces are smooth (the Hausdorff 
dimension for the statistical ensemble of surfaces is two). No rigorous proof of this 
phase transition exists, but numerous numerical results seem to confirm the exis- 
tence of such a transition, and it seems to be a second order transition for the kind 
of extrinsic curvature term we are going to use in this work. In the beginning some 
confusion surrounded this transition. Depending on the details of the discretized 
version of the extrinsic curvature term used, the transition was classified as first, 
second or third order, respectively. It is now understood that only the discretization 
which seems to lead to a second order transition is not "pathological". The other 
actions used, although formally equivalent for smooth surfaces, led to singular sur- 
face configurations if the surface was allowed to fluctuate wildly, as in fact happened 
in the "crumpled phase". This complicated situation highlights the possibility of 
a nontrivial phase structure in random surface theories and the care one has to 
exercise in order to choose a correct discretized action. One can now imagine a 
situation where the attachment of the molecules to specific neighbour molecules in 
the crystalline surface is gradually decreased. The crystalline structure is allowed 
to "melt". Ultimately one will end up with a fluid membrane where the individual 
molecules can move freely. Will the second order phase transition extend all the 
way to a fluid membrane or can one trust the one loop calculation [p| ^[ done for 
the fluid membrane case, which indicated that there is no transition^ ? Computer 
simulations first started by Catterall and later repeated by other groups point to 
the existence of a "crumpling" transition even in the case of fluid membranes, but 
contrary to the situation for crystalline surfaces numerical simulations and analytic 
arguments contradict each other in the fluid membrane case. 

We hope to have convinced the reader that the critical properties of random 
surfaces are a fascinating topic of importance in vastly different areas of physics, 
but that it is an area where only little is known at present. The aim of the present 
article is to develop some theoretical concepts, adequate for the description of the 
critical phenomena of random surfaces, and by extensive numerical simulations try 
to answer some of the questions raised above. 

The paper is organized as follows: In section 2 we discuss the model and define 

3 We thank John Wheater and Frangois David, who independently suggested this interpolating 
scenario to us. 
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the observables, which we use in the canonical numerical simulations. We also relate 
them to the standard grand canonical definitions of the mass gap and string tension. 
Scaling properties of these observables are discussed in section 3. In section 4 we 
describe the system with twisted boundary conditions, which was used in all our 
numerical work to measure the string tension and the mass gap. Numerical results 
are presented in section 5, which is divided into four sub-sections, concerning the 
specific heat, the mass gap, the string tension and the radius of gyration. We 
conclude the paper with a discussion in section 6. 

2 The action and observables 

Let us in this section define the action which we are going to use and discuss the 
observables which will allow us to investigate the critical properties of the system. 
We want to approximate the continuum surfaces discussed in the introduction with 
piecewise flat surfaces. The intrinsic parameter space is then identified with an 
abstract triangulation of appropriate topology, defined entirely by its vertices i and 
links < ij >. The embedding in a D-dimensional euclidean target space is a map 
assigning to each link < ij > a vector X^,fi = 1...D in the target space. The 
numerical results presented in this paper are related solely to the case D = 3. Three 
vectors forming a triangle k should satisfy 



where Xf denotes coordinates of a vertex i in the target space and a sum of vectors 
Xy along any closed path would be zero. For surfaces with handles the embedding 
may be non-trivial and we will discuss this further in section 4. In the rest of this 
section we will assume that the topology is trivial. 

A configuration of the surface is defined by specifying a triangulation (links 
joining the vertices) and by assigning values to vectors X£ satisfying ( |2.1| ). For each 
configuration we define the action S = (3Sg + XSc- Sq is the gaussian part of the 
action 



X% + Xjt k + X£ i = 0. 




If the topology is spherical, fl2.1|) would mean that 



ij j ; 



(2.2) 




<ij> 



(2.3) 



and Sc is the extrinsic curvature part 




(2.4) 



<ij> 
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6ij denoting the angle between the normals to two oriented triangles < ijk > and 
< jil > with a common link < ij >. The coupling constant (3 of the gaussian part 
of the action plays the role of the unit length. 

This action is intended to serve as a discretized version of the terms 

J d 2 ^d a X»d a X» + \J d 2 i^~ d K\i) (2.5) 

in ( p..l|) , where refers to the metric induced by target space. Strictly speaking 
Sc given by (|2.4j ) is not a direct translation of the last term in Q2.5f ), since we are not 
using the induced metric. However, we do not expect there will be any significant 
difference (i.e. the action (|2.5| ) should belong to the correct universality class), since 
both the order of vertices and the shape of triangles in target space are smooth 
functions. 

Two statistical ensembles have our interest: the canonical ensemble where the 
number triangles is kept fixed, and the grand canonical ensemble where the number 
of triangles is allowed to fluctuate. 

The canonical partition function for a closed surface of trivial topology is given 

by 

N 

Zn(P, A) = £ HdX?6(Y,X?) e MS), (2-6) 
TeT i=1 

where the first sum is over all possible triangulations of the surface with N vertices 
and the delta function is necessary because of the translational invariance of the 
action. We can always scale out f3 since Sc is scale invariant and we are left only 
with one coupling constant A. The partition function defines Helmholtz free energy 

F(N,P,\) = -logZ N (/3,\). (2.7) 

Using again the scale invariance of Sc we have the following transformation under 
scaling X? -> yX'?: 

S(X?) = y 2 (3S G (X'^ + XS c {X"t) (2.8) 
and Helmholtz free energy will scale as 

T(N, (3, A) = F(N, y 2 (3, A) - (N - 1)D log y. (2.9) 

The grand canonical ensemble will be defined by 

oo 

Z( l x,P,\)=Y / e-" N Z N (P,\) (2.10) 

N 

where \x is a chemical potential for the number of vertices N (or equivalently the 
number of triangles N2). (3 was not a dynamical coupling constant for the canonical 
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ensemble. The same is true in the grand canonical ensemble where it is clear from 
CT-ffnOD that 

Z(ji,,P,\) = Z(ji l -DlogP,l,\) (2.11) 

It is now obvious that we can always take (3 = 1, and we will do that in the following 
except when explicitly stated differently. The Gibbs free energy can now be defined 
by 

g( fi ,X) = -\ogZ( f i,X) (2.12) 

where we have taken /3 — 1 as announced. 

In the case of the the canonical ensemble we reach the thermodynamic limit by 
taking the size of the system iV — > oo. This thermodynamic limit will in general 
depend on the coupling constant A and for certain values of A there might be phase 
transitions which, in case they are second order, might serve as points where we can 
define a continuum field theory. 

In the case of a grand canonical ensemble we will have a critical line /x = /x c (A) 
in the (A, /x) coupling constant plane, such that the theory is defined for /x > /x c . 
The thermodynamic limit is obtained when number of vertices N diverges, and 
(N) — ► oo corresponds in the language of the grand canonical ensemble to moving 
close to the curve /x = /x c (A) along the /x-axis, starting at large /x. The possibility of 
phase transitions for certain values of A, as described above in the language of the 
canonical ensemble can be addressed in the grand canonical ensemble too. In this 
ensemble such a transition will manifest itself as points on the critical line, which 
separate two different types of critical behaviour when we move along the critical 
line /x = /x c (A). 

For numerical purposes it is much more convenient to work with a canonical 
ensemble. Certain observables are however naturally defined in the grand canonical 
ensemble, and in the following we will show how it is possible to extract information 
about them using only the canonical ensemble. 

The simplest observable which gives us information about the nature of the phase 
transition is the specific heat with respect to A. It can be defined directly in the 
canonical ensemble: 

C(A » " = ^ - <Sc)2 ) (2 ' 13) 

In case there is a second order transition at a finite A c we should see a singularity 
of C(A) as iV — > oo as will be discussed later. 

Even if we observe a second order transition at a finite A c , it does not ensure that 
this transition has anything to do with a string theory. A minimal requirement is 
that the string tension and possible mass excitations scale to zero when A — > A c . The 



following discussion is stimulated by the work of David and Leibler [15]. In order to 
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define the string tension we introduce the grand canonical ensemble of open surfaces 
where the boundary is kept fixed. We imagine the surface will enclose a large area A 
and denote the corresponding partition function Z(fi, A; A). The partition function 
will behave as 

Z(fi, A; A) ~ exp(-(j(yu, X)A), for A -> oo (2.14) 
and we can define the string tension as 

= llm (2.15) 

In the same way we can define the mass gap in the theory as the exponential decay 
of the grand canonical partition function of surfaces with a boundary consisting of 
two points separated by a distance L 

Z(fi, A; L) ~ exp(— m(/i, A)L), for L — > oo (2-16) 

i.e. 

m(Ai,A)= lim G ^ L \ (2.17) 

L—*oo J_j 

In the thermodynamic limit (where /i(A) ~ /•*<;( A) we have (for any boundary 
conditions "-B") the usual relation between Helmholtz free energy T and the Gibbs 
free energy Q 

Q(ji, A; "5") = fiN + J^(AT, A; "5") (2.18) 
iV = ^#^ (2.19) 

If we apply the formula to the case where the boundary condition U B" is the one 
which is used in the definition of the string tension we get 

a(fi, X)A = + JF(iV, A; A) (2.20) 

from which we conclude that 

afa A) = dHN d f A) = v{\ N, A) (2.21) 

with iV expressed in terms \x and A by ( |2.19| ) 

N= da j^ X) A. (2.22) 
on 

In the same way we get where the boundary condition "5" is the one used in the 
definition of the mass gap 

m(/i, \)L = fj,N + F(N, A; L) (2.23) 
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from which we conclude that 



mil*, A) = L) = m(X, N, L), (2.24) 



with N expressed in terms of n and L by ( 12.191) : 



N = dm ^ X \ . (2.25) 



The advantage of \2.2\ ) and j\2. 24 ) is that they allow us to measure the string 



tension and the mass gap using only canonical ensembles. Further, with the action 
in question it follows by a simple scaling argument using (|2.8|) that 

dF(N,\;A) (S G )-±D(N-1) 



dA A 
dF{N,\;L) 2(S G )-D(N 



- (2.26) 

BL L (2 - 27) 

In this way we see that we can measure string tension and mass gap by a simple mea- 
surement of the local observable Sq. To study the critical behaviour it is convenient 
to consider an analogue of the specific heat ( |2.13| ) 

n(XN-A) A " d2 ^N^A) (S G ) — (Sq) 2 — \D(N — 1) 

p(A, N,A) — — — (2.28) 

In deriving (|2.27|) and (|2.28|) we have used only scaling arguments, so the string 
tension should not depend on the shape of the fixed boundary loop. This will only 
be true when the two linear dimensions of the enclosed area are comparable. For 
elongated loops the behaviour will eventually change, since it will be difficult to dis- 
tinguish the physics associated with the area from the (different) physics associated 
with the perimeter. 



3 Scaling properties of the observables 

The aim of the numerical study is to find and understand the critical properties 
of the system. For a canonical system with a fixed number of points N — > oo the 
system is believed to undergo a phase transition for a finite value of the coupling 
constant A. For the target space dimension D = 3 this phase transition seems to be 
of the second order. A standard method to localize the critical point is to study the 
specific heat C(X) (|2.13j ). For a second order phase transition the function C(A) has 
a singularity at A = A c in the thermodynamic limit N ^ oo. For finite values of N, 
C N (\) is a continuous function of A, which develops a maximum at some A = A^. 
Let the maximum of C N (X) be 

CL X = C N (\?). (3.1) 
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From the standard finite-size scaling arguments (e.ff.||30|) we expect for iV — > oo 
that X^ — > \ r with 



|Af-A c | ~ (3.2) 
= A + BN", 

where e and a; are some critical exponents. The increase of C^ ax with growing N is 
considered to be a signal of the second order phase transition unless u = 1, in which 
case the transition could be first order. We will measure C (A) by means of ( |2.13 ). 



However, it is important to stress, as is apparent from ( |2.13|) , that C(A) is not the 
integral of a simple normal-normal correlation function and consequently there is 
no a priori reason that C(X) should diverge even if the normal-normal correlation 
length diverges and in fact our numerical results are compatible with no or less than 
a logarithmic growth of the specific heat with N. 

For the string tension measurement we can use the canonical ensemble of open 
surfaces with a fixed boundary which we will usually take to be a rectangular L\ x L 2 
loop enclosing a given area A, which is not to be confused with the area of the surface. 
As explained in the next section we actually implement the boundary conditions 
differently, but for the discussion of scaling properties the above definition will be 
sufficient. 

The string tension is known not to scale to zero for A = and /i — > /i c (A = 0) 



1% . The thermodynamic limit in the grand canonical ensemble is obtained when 

li R = n - // C (A) -»• (3.3) 
and it is therefore natural to expect a general behaviour 

a(A,/i) = a (A) + d(A)/i^ (A) (3.4) 



with some critical exponent From (p. 22 ), which tells us that N/A = da(fi, X)/dfi 
we then deduce 

N ~ (AM A (3.5) 

It implies that the quantity r = A/N can be expressed in terms of [xr and using eq. 
fl2.2ip we see that the string tension er(A, N, A), defined in the canonical ensemble, 
in the critical region for a fixed A should depend on N and A only through r: 

a(X,N,A)=a(X,r) = ao(X) + d(X)r 2uW/il - 2uW \ (3.6) 



In the derivation of ( |3.6| ) we assumed of course that the finite-size effects can be 
neglected. Formula (|3.6| ) will be used to extract from the numerical simulations the 
critical behaviour of <To(A) and the value of the critical index v(X). The continuum 
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limit is related to the small-r behaviour of ( |3.6| ), where the finite-size corrections 
may be important. This problem will be discussed in the section 5. 

Numerical results on o~(\,r) can be used to get an estimate of the critical ex- 
ponent v, but the behaviour of a(X, p) given by (|3.4j) does not guarantee that the 
continuum theory will have a finite physical string tension. As the thermodynamic 
limit for a fixed A is approached for — > we expect that the possible scaling 
behaviour is extracted from 

a(X, p)A = a phys (X)A phys (3.7) 

where the physical string tension a p h ys (X) and the physical area A phys are kept fixed 
for p^ — > 0. We further want the "bare" area A to go to infinity in order not to deal 
with lattice artifacts. These requirements clearly demand that 00(A) in Q3.4| ) goes 
to zero for A — > A c . Phrased differently the physical string tension will be infinite, 
except at special critical points A c where the coefficient o"o(A c ) vanishes. Assume 
that 

0o(A) ~ (A c - X) a . (3.8) 

(The exponent a can be directly measured in our experiment by fitting the behaviour 
of cr (A) for A — > A c (|3~8| )). We do not expect a scaling of d(\) for A — > A c . In the 
continuum limit we have /ir ; — > and A — > A c with 

A c - A ~ p p R . (3.9) 

The critical exponents a and p satisfy 

ap = 2v, (3.10) 



In ( |3.10|) we assumed that both terms in ( |3.12|) contribute at the same order and we 



now see that the physical area A p h ys in (|3.?1 ) will be related to the "bare" area A by 

A phys ~p%A (3.11) 
and the physical string tension to the "bare" string tension as 

PR PR 

showing explicitly that the physical string tension would be infinite except at the 
critical points where cr (A) vanishes. 

The scaling assignments ( ^.11| ) and fl3.5p are in agreement with the general scaling 



law relation of mass critical exponents like v to the Hausdorff dimension aln of the 
surfaces in target space: 

v=±-. (3.13) 
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In fact one would define the Hausdorff dimension of the present ensemble of surfaces 
by 

N ~ A dH/2 (3.14) 

in the limit where A p h ys is kept fixed but the "bare" A goes to infinity. As we see 
Aphys sets the scale for the divergence of A and iV for {ir ; — > (under the assumption 
that cr (A) = 0): 

A~^, jV~^, (3.15) 



which implies ( |3.13 ): 



N ~ {A) 1/{2u) . (3.16) 

Eq. ( |3.16| ) is not in contradiction to (|3.5| ), but describes how the limit r — > should 
be taken in order to reach the continuum limit. 

As we mentioned in the introduction an additional requirement for an interest- 
ing scaling limit is the correct scaling of other physical observables. An observable 
independent of the string tension is the mass gap, defined in the grand canoni- 
cal ensemble by the exponential decay of the two-point function. In the standard 
approach one fixes two points, separated by a distance y, and sums over all sur- 
faces passing through these points. The transformation of the mass gap definition 
from the grand canonical to the canonical ensemble was discussed in the preceeding 
section (equations fl2.23|) to ( |2.25| )) and we have the relation 

»»> 

Let us now assume that the mass scales at the critical point /x c (A). Opposite 
to the situation for the string tension we expect the bare mass m(A, fx) to scale for 
— fM — /x c (A) — > for all A. Assume the scaling is of the form 

m(A,/i)~c(A)/4 (A) (3.18) 

Using (3.18) one shows that m(X,N,y) defined by ( J2.23 ) to ( |2.25p has a scaling 
behaviour 

m(A, N, y) ~ D(\)t v{X)/{l ~ u{x) \ (3.19) 

where t = y/N is a scaling variable for the mass gap (analogous the r variable for 
the string tension). This formula allows us to determine the critical exponent z/(A) 
from a canonical simulation. 



4 The twisted boundary conditions 

In numerical simulations the systems necessarily have finite dimensions and all mea- 
surements are subject to finite-size effects. In the ideal situation one would like to 
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use these effects to gain additional information about the critical properties of the 
system (like a finite-size scaling analysis). For the string tension measurements one 
may however meet effects of this kind, which are very difficult to estimate and which 
can strongly bias the measured critical behaviour. Following the discussion of the 
previous section, measurements of the string tension would require introducing a 
boundary to the system. For a finite system this may present a serious problem: the 
number of vertices belonging to the boundary may be, and in practice is, a sizeable 
fraction of all vertices. One has also to provide some method of assigning vertices 
to the boundary. No obvious solution seems to be at hand and any solution chosen 
may in fact influence the results. 

In this section we propose a method of avoiding these problems with the help 
of twisted boundary conditions imposed on a surface with the topology of a torus. 
In this case the parameter space of the surface with iV vertices can be visualized as 
a plane, periodic in two non-parallel directions. In each elementary cell we have N 
vertices, each vertex having infinitely many periodic copies. The elementary cells 
can be numbered by two integers k\ and k 2 . Links can connect points both inside 
a single cell and in the neighbouring cells. As before we assign to each link < ij > 
of a lattice in the parameter space a vector X^fi = 1...D, in the target space. We 
assume the map to be periodic. On a torus the embedding may be nontrivial in the 
following sense. As a consequence of ( |2.1|) a sum of vectors along a closed path can 
take values 

E^(n h n 2 ) = n x E^ + n 2 E% (4.1) 

with two constant vectors and E 2 and integers n\ and n 2 denoting the number 
of times the path winds around the two axes of the torus. The fact that the vector 
E^ is non-zero does not contradict the periodicity of X£, because on a torus the 
closed loops with non-zero values of n\ and n 2 can not be contracted to a point. 
We can therefore define non-trivial boundary conditions by choosing two arbitrary 
vectors E% and E 2 . Such a choice of boundary conditions on a torus is the analogue 
of the twisted boundary conditions used in the case of the gauge theories, where 
the gauge fields are periodic up to non-trivial gauge transformations. The vectors 
E£ are topological invariants, they remain unchanged if one or more points along 
the path are shifted in the target space or if the internal geometry of the surface is 
changed by a flip of some of its links. The vector Ej£ does not depend on the point 
index i, but only on the integer distance between the initial and final periods of the 
path. This is a simple consequence of the periodicity of X-i. For nonzero vectors 
Ei and E 2 , ( |4.1| ) means that the coordinates Xf of the vertex % are not strictly 
periodic but depend also on the integer coordinates ki, k 2 numbering the particular 
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elementary cell and are given by 

X?(h, k 2 ) = X? + k x E{ + k 2 E%. (4.2) 

This again does not contradict the periodicity of the dynamical quantities, which 
are functions of X£. 

One can check that this choice of boundary conditions corresponds to spanning 
of the system on a frame E\ x E 2 . This can most easily be seen by considering a 
minimum action configuration of the model. In the embedding space all triangles 
lie in the two-plane spanned by the vectors E% and E 2 with the lengths of the links 
minimizing the Gaussian part of the action. Due to (|4.2|) these lengths are non-zero 
for a non-zero frame. The advantage, as compared to the situation where some 
points are assigned to the boundary is that now the boundary is not present, it 
becomes a translationally invariant concept and the distribution of points is purely 
dynamical. 

In our computation we have chosen the vectors E\ and E 2 to be perpendicular 
with the lengths L\ and L 2 . For the string tension measurements we used the square 
frame L\ = L 2 . For large L\ and L 2 we expect the leading behaviour to be on A = 
LiL 2 . Translational invariance suggests that subleading terms (oc Lx+L 2 ) are absent 
and we should observe only small finite-size corrections of the form Li/L 2 . Another 
possible choice can be L\ ^ and large and L 2 = 0. This clearly corresponds 
to a completely different physical situation, where we measure the analogue of the 
point-point correlation function (or rather loop-loop correlation function) with the 
two loops kept at a distance L\. 

The implementation of the twisted boundary conditions in the numerical simula- 
tions can present some practical problems. In accordance with the discussion above, 
one could be tempted to store the lattice configurations, using the link vectors X& 
rather than the vertex coordinates Xf . Such a parametrization is however danger- 
ous since it allows rounding errors to accumulate. In particular, this can lead to 
a numerical violation of the boundary conditions initially imposed, since both the 
vertex coordinates and the triangulations are changed many times. To avoid this 
problem we decomposed, following ( |4.2|) , the vector X^ as 

X§ = X?-X? + E§, (4.3) 

where = n}-E^ + nf-E? and the integers nf. are non-zero when the link < ij > 
connects points in different elementary cells. To store the configuration we need the 
positions Xj of all points and two integers n}j,nfj for every link. Vectors E%j are 
additive, i.e. + Ejk + Ej-i = for every triangle < ijk >. This implies that also 
the additive. The decomposition (|4.3|) is however not unique. For arbitrary 
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integers l] and If one can perform a transformation 




TP^ L /I TpL 1 I 7 2 Z?^ /I Z?^ / 2 Z?^ 



(4.4) 



which leaves Q4.3| ) invariant 
bounded. 



This transformation can be used to keep the £?y vectors 



In the last section we defined the mass gap by the exponential decay of the two- 
point function, the marked points on the surface separated by a distance y in target 
space. Here we shall consider a different definition, making again use of the twisted 
boundary conditions. This time we consider a system with a zero projected area A p , 
with only one non-zero vector E{. We take E\ = (y, 0, ...), E2 = 0. These boundary 
conditions permit to measure a two-loop function, where two identical loops are 
separated in a translationally invariant way by a distance y and where we sum over 
the loop size. The two-loop function should have the same scaling behaviour as the 
two point function. 

5 Numerical results 

Let us now turn to the numerical simulations of the system. For the measurements 
of the string tension and mass gap we will use the twisted boundary conditions 
described it the last section. For measurements of the specific heat where no ex- 
tended frame in target space is needed these reduce to ordinary periodic boundary 
conditions of surfaces with toroidal topology. We are going to use essentially stan- 
dard Monte Carlo techniques, but before presenting the results, let us collect various 
computational aspects. 

The degrees of freedom in the model are the vertex coordinates and those de- 
scribing the connectivity of the triangulation. Both kinds we update by using a 
local Metropolis algorithm. For the triangulation we use the familiar flip proce- 
dure |Tl|, which is known to be ergodic in the class of all triangulations with the 
same topology. We use the standard restriction that we do not allow non-trivial 
loops of order one and two. In the graph dual to a given triangulation, i.e. a 3 
graph, this restriction corresponds to considering graphs without tadpoles and self- 
energy parts. New values for the coordinates are proposed by performing a heatbath 
or an uj = 2 overrelaxation step J27| with respect to the gaussian part of the ac- 
tion. A Metropolis accept/reject step then takes the curvature term into account. 
We found this coordinate update slightly more efficient than the one in which the 
change in coordinates is generated from a fixed distribution. The gain we measured 
in decorrelation time was somewhat less than a factor two. The flip rate was close 



17 



to 33% in all our calculations. The acceptance rate for the coordinate update varied 
with the simulation parameters and was typically similar to the flip rate. Its lowest 
value was 20%, at A = 1.5 with periodic boundary conditions. 

The dynamical nature of the triangulation restricts the possibilities to vectorize 
the simulation program. We have therefore developed a program in which vector- 
ization is achieved by simulating several systems in parallel. Our results have been 
obtained by simulating 64 parallel systems on a CRAY-YMP. The CPU time re- 
quired per sweep per bond was 5.7 fis. For more details on the algorithm and the 



program, see ref. [26, 32 



While the simultaneous simulation of several systems significantly reduces the 
CPU time needed per sweep, it does, of course, not bring down the number of 
sweeps needed to decorrelate configurations. In fig. 1 we show estimates of the 
integrated autocorrelation time, r int , for the most non-local quantity measured, the 
radius of gyration. n n t refers to the autocorrelation in one system, and the values 
given are averages over the 64 systems considered. We may note that the increase 
in r int between N = 144 and 256, at fixed coupling, corresponds to r int oc N z ' with 
z' roughly 1.6. In our calculations, the number of sweeps used for thermalization 
corresponds in most cases to 20rj nt and always to more than I0r int . Measurements 
were typically taken over 40rj nt . This length of the runs is not very long so one 
might worry about insufficient thermalization. As a check, we therefore carried out 
a few additional much longer runs. The longest one, at (N, A) = (144,1.5) with 
periodic boundary conditions, was two million sweeps, corresponding to 700r int . No 
significant change in the observables was observed in these extended runs, and the 
statistical errors scaled approximately as they should with increasing number of 
iterations. The integrated autocorrelation time did however increase by 15%. The 
numbers given above for r int and its N dependence should therefore be used only as 
a rough guide. 

In the analysis of the specific heat below we use the multi-histogram technique 
by Ferrenberg and Swendsen P§| . This method allows continuation of results ob- 
tained at one or more couplings to, in principle, arbitrary couplings. With limited 
statistics this is, of course, not true, and care is needed in selecting the range for the 
continuation. We made the selection in the following way, similar to that used in 



ref. p8| . We start by determining at each simulated coupling two numbers S\ and 5*2 
such that the probabilities that Sc < Si and that Sc > S 2 are both 25%. We then 
perform a single-histogram continuation, restricted to the coupling interval where 
the resulting value of Sc lies between Si and S 2 . With this restriction we found 
that single-histogram results for the specific heat agreed whenever the couplings 
overlapped. Final numbers were obtained by combining the histograms according 
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to the method in ref. [p9| . The continuation was carried out over a coupling range 
that could be covered by the single-histogram intervals defined earlier. The observed 
consistency between different single-histogram results we take as another indication 
that thermalization effects are under control. 

Finally, we mention that the statistical errors quoted below are jackknife er- 
rors p31 , obtained by taking the results from the different systems as 64 independent 
measurements. 



5.1 Specific heat 

As observed first by Catterall [IB], using spherical topology, the specific heat has a 
maximum at a finite value of the coupling, A ~ 1.5. We have performed a detailed 
study of the specific heat for the case of toroidal topology, with periodic boundary 
conditions. To accurately determine the location and height of the maximum we use 
the multi-histogram technique by Ferrenberg and Swendsen [29|, as described above. 
Fig. 2 shows the specific heat as a function of A for system sizes up to N = 576. 
The curves are the results from multi-histogram analysis, and the points show the 
results at simulated couplings. In agreement with the results from earlier studies 
using spherical topology, we find a maximum in the specific heat near A = 1.5, the 
height of which increases with increasing N. 

The dependence of the position of the maximum, A^ , on the system size is shown 
in fig. 3. The three data points with the largest values of N are well described by 
the form A^ = A c — const N~ a with a = 1/2. However, acceptable fits to this form 
can be obtained for a fairly wide range of a values. As a result of this uncertainty 
about the form of the finite-size corrections, we cannot get a very precise estimate 
of the N —>■ oo limit A c . Fits with x 2 < 1 are obtained for 0.35 < a < 1.2. The 
corresponding bounds on the critical coupling are 1.47 < A c < 1.53. 

How the height of the maximum, C^ ax , varies with the system size is shown in 
fig. 4. The data strongly suggest that C^ ax remains finite in the limit iV — > oo. 
The increase in C^ ax with increasing N becomes slower at large N. The finite-size 
correction to the iV — > oo value appears to vanish faster than 1/N at large N. 



5.2 Mass gap 

In section 2 we saw how one can define a mass gap measurement m(A, N, L) in the 
canonical ensemble. In the thermodynamic limit iV — > oo, this measurement gives 
information about the mass gap m(A, /i), as defined from the exponential decay of 
the grand canonical partition function. We have 

m can (\,t) = lim m(X,N,L) = m(A, //(£)) , (5.1) 

N—>oo,L/N=t 
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where fi(t) is the solution to 

dm(\,n{t)) 1 



(5.2) 



dfi t 

Corresponding to the expected scaling m(A,/i) ~ c(A)/z^ at small we have 
(eq. |3l9| ) m can (\,t) ~ D(A)f (A ) /(1 ~^ A ^ at small t ~ d(A)/4~ KA) . By verifying this 
scaling of m mn (A,t), we can compute the critical exponent z^(A). 

The numerical calculations are restricted to finite N and it is therefore essential 
to keep finite-size effects under control in the estimate of m can (X,t). Let us give an 
estimate of how large iV has to be taken at a given value of t. The finite-size effects 
occur when the asymptotic relation Q(X, /i, L) ~ m(A, /i)L is not fulfilled. This rela- 
tion is valid if terms of lower order in L are negligible, which should be the case when 
m(\, fi)L is large. We thus expect finite-size corrections to be small if m can (\, t)L is 
large. This condition can in the scaling regime be written as t^( A )/( 1_!y ( A ))L ~ 1, or 

tN ~ JV"W (5.3) 

Let us briefly discuss the most likely finite-size effect to the mass gap. The 
generic form of the two point function will be (in the following we will suppress the 
dependence on the coupling constant A) 

Z(iJ,,y)~y-<*e~ m Mv (5.4) 

This means that Gibbs free energy will be of the form 

G{n,y) = m(fi)y + a\ny (5.5) 

and we would expect corrections of the form 

dQ a 

m effW = ~Q^ = m W + ~- ( 5 ' 6 ) 

We can now translate this to the canonical ensemble and get 

m can (t, N) = m can (t) + — . (5.7) 

Our calculations were performed at two different couplings, one in the crumpled 
phase, A = 1.25, and one close to the maximum in the specific heat, A = 1.5. 
Fig. 5 shows the results for m(A, N, L), plotted against t = L/N. In order to check 
the finite-size dependence, we carried out simulations for different N at some fixed 
values of t. As expected from the discussion above, the finite-size corrections turn 
out to be largest at small t. The corrections are large for small N (N = 64 data are 
omitted in the plot), but seem to rapidly become smaller with increasing N. The 
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results for the largest values of N agree well, which suggests that these can be taken 
as reasonable approximations to the thermodynamic limit. 

At large t we find that the results at the two couplings are similar, which is 
expected due to the dominance of the gaussian term in the action. The interesting 
region from the point of view of scaling is, however, at small t, and there we find 
a clear coupling dependence. The results at the three smallest values of t are at 
both couplings consistent with scaling, but with different exponents v(X). Fitting 
the data points with largest N to eq. ( |3.19| ) we get z/(A = 1.25) = 0.279(7) and 



u(X = 1.5) = 0.417(7), with \ 2 near one in both cases. The errors are statistical 
only. Clearly, one would like to confirm that these results reflect the true scaling 
behaviour by going to much smaller values of t. Increasing finite-size effects and 
decorrelation times have prevented us from doing so. In the limited regime probed 
so far, we find that the scaling assumption gives a good description of fairly accurate 
data. 

5.3 String tension 

Our method to calculate the string tension is very similar to that for the mass gap. 
Corresponding to eq. (|5.1|) , we have 



<?can(\r) = lim cr(A, N, A) = <r(A, fi(r)) , (5. 

N— >oo, A/N=r 



where <r(A, N, A) is the measurement defined in section 2 and fi(r) is given by 

<9/i r 



(5.9) 



The major difference is in the expected scaling behaviour. As explained in sec- 
tion 3, we expect <J can (X,r) = cr (X) + a(X)r 2v '( l ~ 2v > with an in general non-zero r 
independent term cr (A). 

In the same way as for the mass gap, we can estimate for which parameter values 
one should expect significant finite-size effects. They should be small if a can (X,r)A 
is large. If we assume the scaling form and that <7o(A) = 0, then we can write this 
condition as r 2v ^ 1 ~ 2v ^A ~ 1, or 

rN ~ N 2u . (5.10) 

Repeating the arguments above we expect the finite-size scaling corrections to be of 
the form 

a can (\, r, N) = a can (X, r) + (5.11) 

rN 

We have carried out simulations at three couplings A = 1.25, 1.4 and 1.5, for 
system sizes up to iV = 576. The results are summarized in fig. 6, where <r(A, N, A) 
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is plotted against r = A/N. Finite-size effects are large for r — > and A — > A c , as 
expected. Therefore, we have not been able to measure the string tension for very 
small r at A = 1.5. It seems, however, that <Jq(X) could vanish near A = 1.5, as 
required for an interesting scaling behaviour. For this we need, in addition, that the 
exponent in 

<WA,r) = <7(A)r w (5.12) 

satisfies uj = 2z/(A)/(l — 2z/(A)), where z/(A) is the mass gap exponent. The value at 
A = 1.5 obtained above, u(X = 1.5) = 0.417(7), corresponds to uj — 4.7 ± 0.5. 

Let us check if the measured string tension at A = 1.5 shows this scaling be- 
haviour. To do so we first extrapolate the results to the thermodynamic limit. In 
fig. 7 we show <r(A = 1.5, A, N) at three fixed values of r for iV = 144,256,400 
and 576. For all three r, we find that the 1/N term gives a statistically acceptable 
description of the finite-size effects, in agreement with ( |5 . 1 1 ) . The lines shown are 
fits to such a form. We can now test for scaling by fitting the extrapolated values to 
eq. Q5.12|) . We find that the data indeed are in good agreement with this behaviour. 
This is illustrated in fig. 8. The exponent obtained, uj = 3.94 ± 0.06, is, moreover, 
close to the value from the mass gap measurement. 

This analysis suggests that there is a coupling close to A = 1.5 where the scaling 
form eq. (|5.12| ) is valid. It is however likely that it is not exactly valid at A = 1.5. We 
therefore want to check how the estimate of uj is affected if we allow for a non-zero 
value of o"o(A = 1.5). To this end we performed a sequence of fits for fixed values of 
uj, having o"o as a free parameter. Fits with y 2 < 1 were obtained for 3.8 < uj < 4.5. 

The determination of uj also involves an assumption about the precise form of 
the finite-size dependence. To get an idea of the importance of this assumption 
we considered the iV = 576 data directly, without any extrapolation procedure. 
Repeating the type of fit with uj fixed and <Jq as a free parameter, we obtained this 
time x 2 < 1 for 3.4 < uj < 4.0. We take this to suggest that the uncertainty in uj 
arising from the N — > oo extrapolation is smaller than 0.5. We then arrive at the 
estimate 3.4 < uj < 5.0, corresponding to 0.38 < u < 0.42. 

The results presented above were all obtained in the crumpled phase or close 
to the phase transition. Unfortunately, we have no conclusive results for the flat 
phase. For a flat surface the minimum in the free energy T is at a finite value of r. 
The canonical string tension is therefore negative for sufficiently small r. We tried 
to verify this by performing some simulations at A = 1.75. Due to large finite-size 
effects and decorrelation times we could however not estimate values for the limit 
of large N. 
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5.4 Radius of gyration 

Another important quantity for the characterization of the surfaces is the radius of 
gyration R g (X,N), which we, in the case of periodic boundary conditions, define as 

Rg(\N) = jfE ( 5 - 13 ) 

i 

The behaviour at large N, 

R g {X,N) ~ N 1 '**™, (5.14) 

is taken to define the Hausdorff dimension d#(A). d#(A) is said to be infinite if the 
growth with N is logarithmic. 

The Hausdorff dimension is expected to decrease with increasing A, from infinity 
at A = to d H = 2 at large A. We here want to estimate the value at the phase 
transition, dff(A c ). To this end we consider the radius of gyration at the maximum 
in the specific heat, R g (X^,N), which we obtain by using again the Ferrenberg- 
Swendsen method. The dependence of R g (X^ , N) on N is shown in the log-log plot 
fig. 9. If we assume that dn(X) varies smoothly across the transition, then the slope 
at large N should give us l/d H (X c ). Over the range of N studied, we find that 
the slope keeps decreasing with increasing N. The two data points with largest 
iV correspond to a value 0^=3.41(7). We take this to suggest the lower bound 
d H (X c ) > 3.4. 



6 Discussion 

We have shown that our data are compatible with a second order transition at 
A = A c ~ 1.5. This transition is characterized by a vanishing string tension a (A) 
for A — > A c and a vanishing mass gap m(A) of the two point function for A — > A c . 
Furthermore 

— ^-4t const. for A — > A c . (6.1) 

m 2 (A) v 1 

This is precisely what we would expect in a string theory without tachyons. It is 
not clear to us what is the effective continuum theory which is compatible with such 
behaviour in three-dimensional target space. As we explained in the introduction 
one part of the motivation for investigating the critical behaviour of strings with 
extrinsic curvature was the relation to fermionic strings where extrinsic curvature 
terms arise in a natural way. From this point of view the results are encouraging: 
Models with extrinsic curvature seem to exhibit non-trivial critical behaviour. From 
the point of view of the superstring it would be interesting to verify that this non- 
trivial behaviour persists for strings embedded in higher dimensions. A further, 
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most interesting problem to address would be the inclusion of the Wess-Zumino-like 
terms mentioned in the introduction. Such terms also arise naturally for fermionic 
strings and could change the critical behaviour even further. It is unfortunately not 
clear to us how these terms can be a part of a numerical study like the one performed 
here since the Wess-Zumino term is complex. 

In the context of conformal field theory the existence of a second order transition 
gives rise to further interesting questions. For the crystalline surfaces we have a well 
defined regular two-dimensional lattice structure, and to each lattice site i we have 
associated d real variables Xi if target space is ci-dimensional. The action between 
these variables is local, although it is not polynomial due to the extrinsic curvature 
term. In case there is a second order transition for a finite value A c of the extrinsic 
curvature coupling one would be tempted to conjecture the existence of an associated 
conformal field theory, characterized by a central charge c. If the extrinsic curvature 
coupling is put to zero we just have a standard gaussian field theory on the lattice 
and it has a gaussian fixed point corresponding to d free fields and a total central 
charge d. In the present context the nature of this fixed point is a little unusual, 
since it is an infrared fixed point and not the usual ultraviolet gaussian fixed point 
known for field theories in less than four dimensions. In fact it seems as if the 
following effective action describes the behaviour of the system for A < A c 



S eff = J d'i x,(0 [-m'(A)^(0 + \d\0\ * M (0 ( 6 - 2 ) 

where 

m(A) =m (A c -A). (6.3) 

Eqs. ( B.2[ ) and (|6.3| ) mean that the long range properties of the system are deter- 
mined by the gaussian part and for this reason one would classify the first fixed 
point away from A = as ultra-violet stable. If we naively apply Zamolodchikov's c- 
theorem we can conclude that the central charge of the conformal theory associated 
with "crumpling" transition has to be larger than or equal d. In fact it is easy to 
calculate the central charge at the fixed point, either analytically or numerically, us- 
ing the effective action ( |6.2|) -( fT3"l) . One gets c = 2d. The numerical results for d = 3 
using the full action for crystalline surfaces, not the approximation indicate 



that the (effective) central charge is between zero and one |fL8|1 , in contradiction with 



the c-theorem. Most likely we are in fact considering a non-unitary conformal field 
theory in which case the c-theorem is formally invalidated. However, one expects a 
c e /j-theorem to apply even in that case and since the central charge measured from 
the partition function in [Tj| is the effective central charge c e ff it does not change the 
apparent contradiction. At least the situation calls for further investigations and one 
should keep in mind that the action involving extrinsic curvature is non-polynomial 
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and it might be that we in this case are mislead by the general principles, which 
do not necessarily apply Whatever the effective central charge is one can couple 
the lattice system to two-dimensional quantum gravity and ask for the properties 
of the corresponding random surface theory. This is what we have done, and since 
the situation already without coupling to gravity (i.e. the crystalline surfaces) ap- 
pears to be non-trivial, the same might be the case when coupled to gravity. In 
fact the connection between the conformal field theory on the world-sheet and the 
scaling properties of the correlation functions in target space has to be non-trivial 
if the model without the Wess-Zumino term will describe a genuine surface theory 
which is euclidean invariant and where the string tension scales. Recent analytic 
results of Polchinski and Strominger indicate that such non-trivial situations can 
indeed occur precisely when the world-sheet variables Xi contain non-polynomial 
terms [35]. Whether our effective string theory with a non-trivial scaling of string 
tension and mass gap has any relation to the scenario suggested in |J5J is presently 
under investigation. 
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Figure Caption 



Fig.l Integrated autocorrelation time for the radius of gyration against system size. 
The boundary conditions are periodic, y\ — y% — 0. 

Fig. 2 The specific heat as a function of A. The bands give the results with errors of 
the multi-histogram analysis. 

Fig.3 Af against iV" 1 / 2 . 

Fig.4 C% ax against 1/N. 

Fig. 5 The mass gap measurement m(A, N, L) plotted against t = L/N. The straight 
lines are fits to eq. ( E>.19| ). 

Fig. 6 The string tension measurement a(X,N,A) plotted against r = A/N. The 
lines connect points with the same A and largest N. 

Fig. 7 <t(A, N, A) against 1/N at three fixed values of r. The lines are fits to the form 
a(X, N, A) = a can (X, r) + c(r)/N. 

Fig. 8 <t(A, N, A) against r u at A = 1.5. The full circles are extrapolations to N — > oo, 
as explained in the text. For the exponent we have used the value u = 3.94, 
obtained from a fit to eq. ( |5.12| ). 

Fig.9 IL(\*,N) against N. 
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